Calcisponges have a ParaHox gene and dynamic expression of dispersed NK homeobox genes

This study was funded by the Sars Centre core budget to M. Adamska. Sequencing was performed at the Norwegian High Throughput Sequencing Centre funded by the Norwegian Research Council. O.M.R. and D.E.K.F. acknowledge support from the BBSRC and the School of Biology, University of St Andrews.Sponges are simple animals with few cell types, but their genomes paradoxically contain a wide variety of developmental transcription factors1,2,3,4, including homeobox genes belonging to the Antennapedia (ANTP) class5,6, which in bilaterians encompass Hox, ParaHox and NK genes. In the genome of the demosponge Amphimedon queenslandica, no Hox or ParaHox genes are present, but NK genes are linked in a tight cluster similar to the NK clusters of bilaterians5. It has been proposed that Hox and ParaHox genes originated from NK cluster genes after divergence of sponges from the lineage leading to cnidarians and bilaterians5,7. On the other hand, synteny analysis lends support to the notion that the absence of Hox and ParaHox genes in Amphimedon is a result of secondary loss (the ghost locus hypothesis)8. Here we analysed complete suites of ANTP-class homeoboxes in two calcareous sponges, Sycon ciliatum and Leucosolenia complicata. Our phylogenetic analyses demonstrate that these calcisponges possess orthologues of bilaterian NK genes (Hex, Hmx and Msx), a varying number of additional NK genes and one ParaHox gene, Cdx. Despite the generation of scaffolds spanning multiple genes, we find no evidence of clustering of Sycon NK genes. All Sycon ANTP-class genes are developmentally expressed, with patterns suggesting their involvement in cell type specification in embryos and adults, metamorphosis and body plan patterning. These results demonstrate that ParaHox genes predate the origin of sponges, thus confirming the ghost locus hypothesis8, and highlight the need to analyse the genomes of multiple sponge lineages to obtain a complete picture of the ancestral composition of the first animal genome.PostprintPeer reviewe

Genome of the ramshorn snail Biomphalaria straminea-an obligate intermediate host of schistosomiasi

This is the final version. Available on open access from Oxford University Press via the DOI in this recordData Availability: The raw genome and RNA sequencing data have been deposited in the SRA under Bioproject No. PRJNA673593. The final chromosome assembly was submitted to NCBI Assembly under accession No. JADKLZ000000000. All data can also be found in the GigaScience Database [75].BACKGROUND: Schistosomiasis, or bilharzia, is a parasitic disease caused by trematode flatworms of the genus Schistosoma. Infection by Schistosoma mansoni in humans results when cercariae emerge into water from freshwater snails in the genus Biomphalaria and seek out and penetrate human skin. The snail Biomphalaria straminea is native to South America and is now also present in Central America and China, and represents a potential vector host for spreading schistosomiasis. To date, genomic information for the genus is restricted to the neotropical species Biomphalaria glabrata. This limits understanding of the biology and management of other schistosomiasis vectors, such as B. straminea. FINDINGS: Using a combination of Illumina short-read, 10X Genomics linked-read, and Hi-C sequencing data, our 1.005 Gb B. straminea genome assembly is of high contiguity, with a scaffold N50 of 25.3 Mb. Transcriptomes from adults were also obtained. Developmental homeobox genes, hormonal genes, and stress-response genes were identified, and repeat content was annotated (40.68% of genomic content). Comparisons with other mollusc genomes (including Gastropoda, Bivalvia, and Cephalopoda) revealed syntenic conservation, patterns of homeobox gene linkage indicative of evolutionary changes to gene clusters, expansion of heat shock protein genes, and the presence of sesquiterpenoid and cholesterol metabolic pathway genes in Gastropoda. In addition, hormone treatment together with RT-qPCR assay reveal a sesquiterpenoid hormone responsive system in B. straminea, illustrating that this renowned insect hormonal system is also present in the lophotrochozoan lineage. CONCLUSION: This study provides the first genome assembly for the snail B. straminea and offers an unprecedented opportunity to address a variety of phenomena related to snail vectors of schistosomiasis, as well as evolutionary and genomics questions related to molluscs more widely

Evolution through segmental duplications and losses : A Super-Reconciliation approach

The classical gene and species tree reconciliation, used to infer the history of gene gain and loss explaining the evolution of gene families, assumes an independent evolution for each family. While this assumption is reasonable for genes that are far apart in the genome, it is not appropriate for genes grouped into syntenic blocks, which are more plausibly the result of a concerted evolution. Here, we introduce the Super-Reconciliation problem which consists in inferring a history of segmental duplication and loss events (involving a set of neighboring genes) leading to a set of present-day syntenies from a single ancestral one. In other words, we extend the traditional Duplication-Loss reconciliation problem of a single gene tree, to a set of trees, accounting for segmental duplications and losses. Existency of a Super-Reconciliation depends on individual gene tree consistency. In addition, ignoring rearrangements implies that existency also depends on gene order consistency. We first show that the problem of reconstructing a most parsimonious Super-Reconciliation, if any, is NP-hard and give an exact exponential-time algorithm to solve it. Alternatively, we show that accounting for rearrangements in the evolutionary model, but still only minimizing segmental duplication and loss events, leads to an exact polynomial-time algorithm. We finally assess time efficiency of the former exponential time algorithm for the Duplication-Loss model on simulated datasets, and give a proof of concept on the opioid receptor genes

Are Hox Genes Ancestrally Involved in Axial Patterning? Evidence from the Hydrozoan Clytia hemisphaerica (Cnidaria)

Background: The early evolution and diversification of Hox-related genes in eumetazoans has been the subject of conflicting hypotheses concerning the evolutionary conservation of their role in axial patterning and the pre-bilaterian origin of the Hox and ParaHox clusters. The diversification of Hox/ParaHox genes clearly predates the origin of bilaterians. However, the existence of a "Hox code' predating the cnidarian-bilaterian ancestor and supporting the deep homology of axes is more controversial. This assumption was mainly based on the interpretation of Hox expression data from the sea anemone, but growing evidence from other cnidarian taxa puts into question this hypothesis. Methodology/Principal Findings: Hox, ParaHox and Hox-related genes have been investigated here by phylogenetic analysis and in situ hybridisation in Clytia hemisphaerica, an hydrozoan species with medusa and polyp stages alternating in the life cycle. Our phylogenetic analyses do not support an origin of ParaHox and Hox genes by duplication of an ancestral ProtoHox cluster, and reveal a diversification of the cnidarian HOX9-14 genes into three groups called A, B, C. Among the 7 examined genes, only those belonging to the HOX9-14 and the CDX groups exhibit a restricted expression along the oralaboral axis during development and in the planula larva, while the others are expressed in very specialised areas at the medusa stage. Conclusions/Significance: Cross species comparison reveals a strong variability of gene expression along the oral-aboral axis and during the life cycle among cnidarian lineages. The most parsimonious interpretation is that the Hox code, collinearity and conservative role along the antero-posterior axis are bilaterian innovations

SynBlast: Assisting the analysis of conserved synteny information

<p>Abstract</p> <p>Motivation</p> <p>In the last years more than 20 vertebrate genomes have been sequenced, and the rate at which genomic DNA information becomes available is rapidly accelerating. Gene duplication and gene loss events inherently limit the accuracy of orthology detection based on sequence similarity alone. Fully automated methods for orthology annotation do exist but often fail to identify individual members in cases of large gene families, or to distinguish missing data from traceable gene losses. This situation can be improved in many cases by including conserved synteny information.</p> <p>Results</p> <p>Here we present the <monospace>SynBlast</monospace> pipeline that is designed to construct and evaluate local synteny information. <monospace>SynBlast</monospace> uses the genomic region around a focal reference gene to retrieve candidates for homologous regions from a collection of target genomes and ranks them in accord with the available evidence for homology. The pipeline is intended as a tool to aid high quality manual annotation in particular in those cases where automatic procedures fail. We demonstrate how <monospace>SynBlast</monospace> is applied to retrieving orthologous and paralogous clusters using the vertebrate <it>Hox </it>and <it>ParaHox </it>clusters as examples.</p> <p>Software</p> <p>The <monospace>SynBlast</monospace> package written in <monospace>Perl</monospace> is available under the GNU General Public License at <url>http://www.bioinf.uni-leipzig.de/Software/SynBlast/</url>.</p

Coordinated spatial and temporal expression of Hox genes during embryogenesis in the acoel Convolutriloba longifissura

Background: Hox genes are critical for patterning the bilaterian anterior-posterior axis. The evolution of their clustered genomic arrangement and ancestral function has been debated since their discovery. As acoels appear to represent the sister group to the remaining Bilateria (Nephrozoa), investigating Hox gene expression will provide an insight into the ancestral features of the Hox genes in metazoan evolution. Results: We describe the expression of anterior, central and posterior class Hox genes and the ParaHox ortholog Cdx in the acoel Convolutriloba longifissura. Expression of all three Hox genes begins contemporaneously after gastrulation and then resolves into staggered domains along the anterior-posterior axis, suggesting that the spatial coordination of Hox gene expression was present in the bilaterian ancestor. After early surface ectodermal expression, the anterior and central class genes are expressed in small domains of putative neural precursor cells co-expressing ClSoxB1, suggesting an evolutionary early function of Hox genes in patterning parts of the nervous system. In contrast, the expression of the posterior Hox gene is found in all three germ layers in a much broader posterior region of the embryo. Conclusion: Our results suggest that the ancestral set of Hox genes was involved in the anteriorposterior patterning of the nervous system of the last common bilaterian ancestor and were later co-opted for patterning in diverse tissues in the bilaterian radiation. The lack of temporal colinearity of Hox expression in acoels may be due to a loss of genomic clustering in this clade or, alternatively, temporal colinearity may have arisen in conjunction with the expansion of the Hox cluster in the Nephrozoa

Changing Hydrozoan Bauplans by Silencing Hox-Like Genes

Regulatory genes of the Antp class have been a major factor for the invention and radiation of animal bauplans. One of the most diverse animal phyla are the Cnidaria, which are close to the root of metazoan life and which often appear in two distinct generations and a remarkable variety of body forms. Hox-like genes have been known to be involved in axial patterning in the Cnidaria and have been suspected to play roles in the genetic control of many of the observed bauplan changes. Unfortunately RNAi mediated gene silencing studies have not been satisfactory for marine invertebrate organisms thus far. No direct evidence supporting Hox-like gene induced bauplan changes in cnidarians have been documented as of yet. Herein, we report a protocol for RNAi transfection of marine invertebrates and demonstrate that knock downs of Hox-like genes in Cnidaria create substantial bauplan alterations, including the formation of multiple oral poles (“heads”) by Cnox-2 and Cnox-3 inhibition, deformation of the main body axis by Cnox-5 inhibition and duplication of tentacles by Cnox-1 inhibition. All phenotypes observed in the course of the RNAi studies were identical to those obtained by morpholino antisense oligo experiments and are reminiscent of macroevolutionary bauplan changes. The reported protocol will allow routine RNAi studies in marine invertebrates to be established

An Amphioxus Gli Gene Reveals Conservation of Midline Patterning and the Evolution of Hedgehog Signalling Diversity in Chordates

Background. Hedgehog signalling, interpreted in receiving cells by Gli transcription factors, plays a central role in the development of vertebrate and Drosphila embryos. Many aspects of the signalling pathway are conserved between these lineages, however vertebrates have diverged in at least one key aspect: they have evolved multiple Gli genes encoding functionally-distinct proteins, increasing the complexity of the hedgehog-dependent transcriptional response. Amphioxus is one of the closest living relatives of the vertebrates, having split from the vertebrate lineage prior to the widespread gene duplication prominent in early vertebrate evolution. Principal findings. We show that amphioxus has a single Gli gene, which is deployed in tissues adjacent to sources of hedgehog signalling derived from the midline and anterior endoderm. This shows the duplication and divergence of the Gli family, and hence the origin of vertebrate Gli functional diversity, was specific to the vertebrate lineage. However we also show that the single amphioxus Gli gene produces two distinct transcripts encoding different proteins. We utilise three tests of Gli function to examine the transcription regulatory capacities of these different proteins, demonstrating one has activating activity similar to Gli2, while the other acts as a weak repressor, similar to Gli3. Conclusions. These data show that the vertebrates and amphioxus have evolved functionally-similar repertoires of Gli proteins using parallel molecular routes; vertebrates via gene duplication and divergence, and amphioxus via alternate splicing of a single gene. Our results demonstrate that similar functional complexity of intercellular signalling can be achieved via different evolutionary pathways

Hox cluster duplication in the basal teleost Hiodon alosoides (Osteoglossomorpha)

Large-scale—even genome-wide—duplications have repeatedly been invoked as an explanation for major radiations. Teleosts, the most species-rich vertebrate clade, underwent a “fish-specific genome duplication” (FSGD) that is shared by most ray-finned fish lineages. We investigate here the Hox complement of the goldeye (Hiodon alosoides), a representative of Osteoglossomorpha, the most basal teleostean clade. An extensive PCR survey reveals that goldeye has at least eight Hox clusters, indicating a duplicated genome compared to basal actinopterygians. The possession of duplicated Hox clusters is uncoupled to species richness. The Hox system of the goldeye is substantially different from that of other teleost lineages, having retained several duplicates of Hox genes for which crown teleosts have lost at least one copy. A detailed analysis of the PCR fragments as well as full length sequences of two HoxA13 paralogs, and HoxA10 and HoxC4 genes places the duplication event close in time to the divergence of Osteoglossomorpha and crown teleosts. The data are consistent with—but do not conclusively prove—that Osteoglossomorpha shares the FSGD

Classification and nomenclature of all human homeobox genes

<p>Abstract</p> <p>Background</p> <p>The homeobox genes are a large and diverse group of genes, many of which play important roles in the embryonic development of animals. Increasingly, homeobox genes are being compared between genomes in an attempt to understand the evolution of animal development. Despite their importance, the full diversity of human homeobox genes has not previously been described.</p> <p>Results</p> <p>We have identified all homeobox genes and pseudogenes in the euchromatic regions of the human genome, finding many unannotated, incorrectly annotated, unnamed, misnamed or misclassified genes and pseudogenes. We describe 300 human homeobox loci, which we divide into 235 probable functional genes and 65 probable pseudogenes. These totals include 3 genes with partial homeoboxes and 13 pseudogenes that lack homeoboxes but are clearly derived from homeobox genes. These figures exclude the repetitive <it>DUX1 </it>to <it>DUX5 </it>homeobox sequences of which we identified 35 probable pseudogenes, with many more expected in heterochromatic regions. Nomenclature is established for approximately 40 formerly unnamed loci, reflecting their evolutionary relationships to other loci in human and other species, and nomenclature revisions are proposed for around 30 other loci. We use a classification that recognizes 11 homeobox gene 'classes' subdivided into 102 homeobox gene 'families'.</p> <p>Conclusion</p> <p>We have conducted a comprehensive survey of homeobox genes and pseudogenes in the human genome, described many new loci, and revised the classification and nomenclature of homeobox genes. The classification scheme may be widely applicable to homeobox genes in other animal genomes and will facilitate comparative genomics of this important gene superclass.</p